Breathing modes of a fast rotating Fermi gas 
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We derive the frequency spectrum of the lowest compressional oscillations of a 3D harmonically 
trapped Fermi superfluid in the presence of a vortex lattice, treated in the diffused vorticity approx- 
imation within a hydrodynamic approach. We consider the general case of a superfluid at T = 
characterized by a polytropic equation of state (~ n 7 ), which includes both the Bose-Einstein con- 
densed regime of dimers (7 = 1) and the unitary limit of infinite scattering length (7 = 2/3). 
Important limiting cases are considered, including the centrifugal limit, the isotropic trapping and 
the cigar geometry. The conditions required to enter the lowest Landau level and quantum Hall 
regimes at unitarity are also discussed. 
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The experimental realization of quantized vortices in 
interacting ultracold Fermi gases |l| has opened new chal- 
lenging perspectives in the experimental and theoretical 
study of superfluidity. These perspectives are particu- 
larly important because in Fermi gases the measurement 
of the order parameter is not directly accessible. Actu- 
ally, only for small and positive values of the scattering 
length, when dimers built up with pairs of atoms of op- 
posite spin are formed, the fermionic system gives rise to 
the phenomenon of Bose-Einstein condensation (BEC), 
whose onset is clearly revealed by the bimodal structure 
of the density Q . The observation of vortices for negative 
values of the scattering length as well as in the unitary 
limit close to a Feshbach resonance, where the scattering 
length is larger than the interparticle distance, conse- 
quently provides a unique source of information on the 
superfluid nature of these novel configurations. 

The measurements of Ref. [l| have however shown that 
vortices in atomic Fermi gases are not directly observable 
in situ nor after expansion, unless one suddenly ramps 
the scattering length to small and positive values (cor- 
responding to the BEC regime) just after the release of 
the trap. Indeed, the visibility of vortices is limited in 
both the BCS and unitary limits, in the first case due to 
the reduced contrast and in the latter mainly due to the 
smallness of their size, fixed by the interparticle distance. 

For the above reasons it is interesting to explore more 
macroscopic signatures of the presence of vortices. A first 
important source of information comes from the bulge ef- 
fect associated with the increase of the radial size of the 
cloud produced by the centrifugal force. More system- 
atic information comes from the study of the collective 
oscillations. For example the splitting of the quadrupole 
frequencies with opposite angular momentum [3(] provides 
direct information on the angular momentum carried by 
the vortical configuration and has been used to measure 
even the quantization of a single vortex line in Bose- 
Einstein condensed atomic gases Q. In this work we 



focus on the study of the compressional modes whose 
frequency is affected by the presence of the vortex lines 
and, at the same time, is sensitive to the nature of the 
configuration (BEC gas of dimers, unitary limit, etc.). 
The study of the compressional modes actually provides 
a unique information on the equation of state of these 
systems. In the absence of rotation it has been the ob- 
ject of recent theoretical d, Q and experimental @, H, @| 
work. Theoretical studies of the collective oscillations in 
rotating Fermi gases were so far limited to small angular 
velocities in the absence of vortices [lfjj. First calcula- 
tions based on the diffused vorticity approximation were 
recently carried out in cylindrical geometry 11]. 

We consider here a two-component Fermi gas with bal- 
anced spin population, trapped by a three-dimensional 
axisymmetric harmonic potential, so that the collective 
oscillations can be labeled by the axial component m of 
angular momentum. The typical wavelength of the low- 
est modes is of the order of the system size. When the 
number of vortices in the sample is large, this length 
scale is much larger than the intervortex distance. In 
order to evaluate the corresponding oscillation frequen- 
cies one can then rely on a coarse grain description of 
the syjrtern, the so-called diffused vorticity approxima- 
tion [l2, [l3(, which does not require to deal with the 
microscopic details of single vortices. By considering the 
case of a regular lattice of singly quantized vortices, this 
long-wavelength description assumes the vorticity to be 
uniformly spread in the fluid. In practice, if the vortex 
lattice rotates at angular velocity H, the average curl of 
the velocity field v is given by V Ao = 2fi, characterizing 
the rigid body rotation v — flAr. This also corresponds 
to a uniform vortex density, which, for a Fermi super- 
fluid, is n v = 2MD,/%H, i.e., a factor 2 larger than in 
the case of a Bose superfluid with the same value of the 
atomic mass M. 

The main consequence of the diffused vorticity approx- 
imation is the introduction of an effective velocity field 
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which does not satisfy anymore the irrotationality con- 
straint of the microscopic superfluid flow, but accounts 
for the presence of the vortex lattice. Within this frame- 
work, we consider the problem of solving the equations 
of rotational hydrodynamics with a polytropic equation 
of state, where the chemical potential fi is assumed to 
have a power law dependence on the density n, namely, 
fi oc n 1 . This parametrization treats exactly several 
important configurations of interacting Fermi gases, in- 
cluding the Bose-Einstein condensed regime of dimcrs 
(7 = 1), where the scattering length is small and pos- 
itive, and the unitary regime of infinite scattering length 
where the equation of state takes a universal density de- 
pendence characterized by the value 7 = 2/3. Moreover, 
as discussed in the literature (see, e.g., Ref. Q and ref- 
erences therein), it is possible to provide an accurate, al- 
though approximate, description of the entire BEC-BCS 
crossover by introducing an effective exponent 7 for the 
equation of state. Since the compressional modes are 
sensitive to the equation of state, the accurate study of 
their frequency can then provide a useful insight on the 
various regimes achieved in the experiments. 

The equations of rotational hydrodynamics, written, 
in the laboratory frame, are given by 



d t n + V • (nv) = , 



(1) 



Md t v + V (Mv 2 /2 + V cxt + /iioc) = Mv A (V A v) ,(2) 

where u\ oc (r,t) oc n 7 (r,i) is the local chemical poten- 
tial fixed by the equation of state of the uniform matter, 
v(r,t) is the velocity field, and Kxt is the external po- 
tential which is assumed to be the same for both the spin 
components of the Fermi gas. For an axisymmetric har- 
monic potential V cyLt = M \ui"]_{x 2 + y 2 ) + lj 2 z 2 ] /2 and 
for a rotation of the trap in the x-y plane at frequency 
Slo, the equilibrium solutions of Eqs. (HJ and |J2J) are given 
by v (r) = O A r and n (r) oc [fj, - V c ^_t{r)] lh , where 
Kxt = Vexx — MQq(x 2 + y 2 )/2 is the renormalized trap- 
ping potential accounting for the centrifugal effect pro- 
duced by the rotation and the chemical potential uq is 
obtained from the normalization condition for the den- 
sity. The centrifugal force causes a bulge effect which 
modifies the aspect ratio of the rotating cloud according 
to the relationship 

r 2 jr\ = (wi - ng) A4 (3) 

where R z and R± are, respectively, the axial and radial 



Thomas-Fermi radii of the cloud. It also fixes a natural 
limit for the angular velocity f2o which cannot exceed the 
radial trapping frequency lu±. 

By expanding Eqs. |T]) and (J2J) with respect to small 
perturbations of the density and velocity field, n — 
no + Sn and v = Vq + Sv, one obtains two coupled lin- 
earized equations which admit several solutions of rele- 
vant physical interest. On the one side, one has surface 
solutions carrying angular momentum and characterized 
by irrotational flow. For example, the most relevant 
m = ±2 quadrupole solutions are described by density 
variations of the form Sn oc (x ± ry) 2 and by the veloc- 
ity field 5vocV(x± iy) 2 [U EDI, B Ej- Thes e surface 
modes exhibit the dispersion ui± = y // 2uj 2 _ — 0§ ± f2o and 
are strongly affected by the rotational effect, as experi- 
mentally proven in the case of Bose-Einstein condensed 
atomic gases [l6j |. Their frequency is however indepen- 
dent of the equation of state and is not expected to ex- 
hibit a new behavior in the case of a Fermi superfluid. 

In addition to the surface modes the hydrodynamic 
equations exhibit an important class of m = compres- 
sional modes. In order to solve the linearized equations 
of motion we use the Ansatz 1(J HI 



Sv = {SO Ar + V [a ± {x 2 + y 2 ) + a z z 2 ]} e~ lujt ,(4) 



Sn = rip 7 loo + a±(x 2 + y z ) + a z z z | t 



(5) 



where <5J7, parallel to the axial direction, accounts for 
the proper variation of the angular velocity during the 
oscillation. The inclusion of this term is crucial to ensure 
the conservation of an gula r momentum. It was ignored 
in previous works llj, LLJ] on the collective oscillations 
of two-dimensional Bose and Fermi gases containing a 
vortex lattice where pure irrotational flow was assumed. 
While the irrotational assumption is valid for the m/0 
modes, it turns out to be inadequate for the compres- 
sional m = oscillations which are associated with vari- 
ations of the vortex density and hence of the angular 
velocity. 

The Ansatz ([4j|5]) gives rise to a linear system yielding 
three solutions for the oscillation frequency ui. One of 
them is the trivial solution lu = 0, corresponding to a 
change of the equilibrium configuration due to the adia- 
batic change of the angular velocity of the system. The 
other two solutions instead correspond to the radial and 
axial breathing modes of the gas and their frequency is 
given by 



(l + 7)wi + 



2 + 7 



lo 2 z + (1 - 7 )ng 



± 



^(1 + 7 )M + (^P) 2 "i + (1 - 7) 2 ^ + (7 2 - 3 7 - 2)^K + 2(1 + 7 )(1 - 7Kfi§ - (7 2 ~ 7 + 2) W ^g . (6) 
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FIG. 1: Frequencies of the breathing modes lu-^/lu± (a) and 
u)—/uj±_ (b) of Eq.j6}, as a function of the rotation frequency 
Qo/uj±. We consider a cigar-shape trapping (u>j_/ui z = 10), and 
show results for the BEC (7 = 1) and the Unitary (7 = 2/3) 
regimes. 



Equation ^ represents the main result of the present 
work. These two solutions arise from the coupling be- 
tween the radial and axial motion caused by the hydro- 
dynamic forces in the presence of the rigid rotation of the 
gas. It is easy to see that Eq. ([6]) reproduces, as limiting 
cases, the result of Ref. [lTj for the non-rotating Bose 
condensed gas (fin = 0, 7 = 1), the result of Ref. [l3| for 
the rotational Bose gas (On 7^ 0, 7 = 1), and finally the 
result of Ref. [10| for the non-rotating superfluid Fermi 
gas (On = 0, for generic 7). 

The predicted behavior of the frequencies for a typical 
cigar-shaped geometry (uj±/ui z = 10) is shown in Fig[TJ 
where we explicitly compare the BEC (7 = 1) and the 
unitary (7 = 2/3) regimes. 

Let us now discuss some limiting cases predicted by 
Eq. (J6j> . A first important case is given by the centrifugal 
limit Slo — > wj_ where the cloud assumes a disk shape as a 
consequence of the bulge effect ([3]) and the two solutions 



<j6j) take the simple form 

id 1 = 2 id I 



v/2 + 7< 



(7) 
(8) 



In the centrifugal limit the frequency of the radial breath- 
ing mode approaches the universal value 2td±, indepen- 
dent of the equation of state [see Fig- dU)] and of the 
value of trap deformation lu±/lu z , while the 7 depen- 
dent frequency of the axial breathing mode coincides with 
the value predicted by the hydrodynamic equations at 
On = 0, in the disk-shaped configuration lo±_ <C oj z - In 
this regard, one should notice that result JHJ), as well as 
the more general result ([6]), has been derived assuming a 
3D configuration, i.e., assuming the validity of the local 
density approximation along the three directions. When 
On becomes too close to u± the gas becomes extremely di- 
lute and the Thomas-Fermi condition /in 3> hcd z is even- 
tually violated with the consequent transition to a 2D 
configuration. In this case the axial frequency takes the 
ideal gas value 2id z instead of %J2 + 7 uj z . For a Bose- 
Einstein condensed atomic gas (7 = 1) this transition 
has been investigated experimentally |18| . In the case of 
a Fermi gas at unitarity the conditions required to reach 
the 2D regime by approaching the centrifugal limit are 
much more severe (see discussion below) . 

Another interesting configuration is given by 



isotropic trap geometry ui z 
reduces to 



for which Eq. 



the 
© 



•37 



+ (l- 7 ) ft 



W~7 2 



(1 



7) 2 



7(37-1)^ . (9) 



It is worth noticing that, while in the absence of rota- 
tion the frequency w_ reduces to the result V^coq for the 
surface quadrupole m = mode and the frequency lo + 
approaches the value \J2 + 3jujq of the pure monopole 
compression mode, the rotation provides a coupling be- 
tween the two modes even for isotropic trapping, so that 
they are both sensitive to the value of the equation of 
state. A special case is the unitary regime (7 = 2/3) 
where the two solutions reduce to 



= 2 LdQ , 



2^o 



()'-! 

"0 



(10) 
(11) 



The spherical trapping geometry in the unitary regime 
is actually of particular interest since in this case the 
Schrodinger equation exhibits important scaling proper- 
ties [l9| . These give rise to universal features for the free 
expansion as well as for the radial monopole frequency 
lu + , which turns out to be independent of f2o- It is worth 
stressing that, because of the centrifugal effect, the shape 
of the gas is not spherical in spite of the spherical symme- 
try of the trap. The dynamics are however isotropic, the 
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solution of the equations of motion being exactly fixed 
by an isotropic scaling transformation. 

Let us now study the experimentally relevant case of a 
strongly anisotropic trap tu z <C (cigar shape) . In this 
case Eq. ^ yields the useful results [20| 



u+ = ^2(1 + 7) (4 + 2(1-7) fig , 

/ 2 + 3 7 +(2- 7 ) (r> /^ ± ) 2 
^ V l + 7 +(l-7)(«o/^) 2 ' 



(12) 
(13) 



which correspond to the solution for the radial and ax- 
ial breathing modes, respectively. It is remarkable to see 
that in the BEC case (7 = 1) the frequency of the radial 
breathing mode is independent of fio, reflecting the pecu- 
liar behavior exhibited by the Gross-Pitaevskii equation 
in 2D 21]. Conversely, in the unitary regime (7 = 2/3), 
the previous equations reduce to 



(14) 
(15) 




(fio/fJ-)' 



showing that the radial breathing mode depends on the 
actual value of the angular velocity and its frequency 
ranges from the value y/W/3cj± at Qq = (recently mea- 
sured in Lithium gases @; B @j) to the universal value 
2lo± holding in the centrifugal limit [see Fig.((T])]. 

It is useful to recall that the hydrodynamic description 
employed in the present work is based on the applicability 
of the Thomas-Fermi approximation fiQ 3> huj±_ which, in 
the presence of vortices implies that the vortex size, fixed 
by the healing length £, be much smaller than the inter- 
vortex distance d v = y/h/Mfl: 



(16) 



When £ ~ d v one enters the lowest Landau level (LLL) 
regime, a regime already explored in rotating Bose gases 
both theoretically [23| and experimentally [l8| . The con- 
dition no ^> hui, instead ensures the 3D nature of the 
configuration [22j . At unitarity, where the healing length 
is of the order of the interparticle distance d, the transi- 
tion to the LLL regime (£ ~ d v ) is reached for angular 
velocities extremely close to the centrifugal limit, satis- 
fying the condition 



1 

N' 



(17) 



Here, for simplicity, we have considered isotropic trap- 
ping (u>± ~ lo z = wo), so that the condition (fT7|) also 
implies the transition to 2D. It is immediate to see that 
in a Fermi gas at unitarity this condition is equivalent 
to requiring that the number of vortices be close to the 



number of particles (N v ~ N) [24[ , a regime where one 
expects to observe quantum Hall (QH) like effects. 

It is worth comparing Eq. (117p with the condition re- 
quired to reach the LLL and QH regimes in a dilute Bose 
gas. In this case the LLL regime (also corresponding to 
the 2D regime) is obtained under the less severe require- 
ment [1 — (SI/ujo) 2 ] — 1/G, where G — Na/aho- Indeed, 
since the scattering length a is much smaller than the os- 
cillator length a/i = y/h/Mu)Q one has G <C N. The con- 
dition for reaching the quantum Hall regime (N v ~ N) 
is instead given by the much more severe requirement 

[1- (n /lo ) 2 } ^g/n\ 

An other interesting case concerns the BCS regime 
of negative values of the scattering length. As recently 
pointed out in Rcf. 25], in this case the centrifugal limit 
cannot be reached by keeping the system in the super- 
fluid phase. In fact, since the density of the gas becomes 
smaller and smaller as il — > lo±, the pairing gap even- 
tually becomes of the order of the trapping frequency 
and superfluidity is lost, with the emergence of a smooth 
transition between a supcrfluid central core containing a 
vortex lattice and a rotating normal fluid at the periph- 
ery While the equations of rotational hydrodynamics are 
expected to hold also in the normal phase, the detailed 
structure of elementary excitations might be influenced 
by the co-existence of the normal and superfluid compo- 
nents. 

Let us finally point out that the equations of hydrody- 
namics are well suited also to study the problem of the 
expansion and in particular to predict how the presence 
of vorticity affects the time evolution of the aspect ratio 
after release of the trap. 

In conclusion we have derived the frequency spectrum 
of the breathing modes of a 3D harmonically trapped 
Fermi superfluid in the presence of a vortex lattice. Spe- 
cial attention has been devoted to the unitarity regime, 
where the collective frequencies are found to exhibit a dif- 
ferent dependence on the angular velocity with respect to 
the case of a dilute Bose gas. 
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